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Estimation of stellar parameters using Monte 
Carlo Markov Chains. 
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Abstract. We apply Monte Carlo Markov Chain methods to the stellar parameter estima- 
tion problem. This technique is useful when dealing with non-linear models and allows to 
derive realistic error bars on the inferred parameters. We give the first results obtained for 
a Cen A. 
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1. Introduction 

Estimation of the parameters of stellar mod- 
els is a challenging task. The increasing pre- 
cision achieved in the determination of atmo- 
spheric parameters, coupled to the fast growing 
amount of available seismic and interferomet- 
ric data, has allowed to improve the constraints 
on stellar models. 

However, it is also important to improve 
the methods for stellar parameter estimation. 
Stellar models are known to be highly non- 
linear. Therefore, optimisation methods may 
fail in computing the optimum parameters. 
Moreover, it is usually difficult to associate a 
confidence level to the inferred stellar parame- 
ters. 

We consider here the use of Monte Carlo 
Markov Chain (MCMC) algorithms to address 
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this problem. They offer the advantage of be- 
ing efficient for non-linear models and allow, 
when combined with statistical inference, to 
estimate jointly the parameters and associated 
confidence levels. 



2. Principle of the MCMC algorithm 

Let 9 and X be the vectors collecting respec- 
tively the model parameters and the observa- 
tional data. Our probabilistic approach aims 
at estimating the posterior probability distri- 
bution (PPD) of the model parameters, condi- 
tional on the available data for the star, which 
is given by Bayes's formula : 



7r{e\X) = 



K ' 



(1) 
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with 7r(6i|X) the PPD, /(X|0) the likelihood, 
7r{9) a given prior density function for the 
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Table 1. Results of the MCMC simulations for a Cen A. For each run the first row gives the mean 
value and the corresponding 1-cr error bar. The second row gives the parameter value estimated 
from the maximum a posteriori (MAP, in brackets). 



Run* 


M(Mo) 


Tm (Mo) 


t. (Gyr) 


o-u (Gyr) 


1 


1.122 


0.019 


5.59 


1.13 




(1.119) 




(5.80) 




2 


1.107 


0.007 


6.52 


0.45 




(1.106) 




(6.55) 




3 


1.117 


0.014 


5.91 


0.51 




(1.116) 




(5.97) 





model parameters and K a normalisation con- 
stant. Assuming that the observations are inde- 
pendent, we used a Gaussian likelihood: 



f{X\e) cc exp 



1 ^ 
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xl\0) - Xi 



o-i 
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with Xf^{6) the output of the model correspond- 
ing to the /-th component of X, cr,- the standard 
deviation, chosen as the Icr-error bar on the 
measurement and the number of measure- 
ments. 

The idea behind MCMC algorithms is to 
generate samples distributed according to a 
distribution of interest, here n{0\X). Since these 
distributions are usually complex, samples are 
generated using simpler ones, called instru- 
mental distributions. The useful property of an 
MCMC algorithm is that when it reaches con- 
vergence it generates sets of parameters ac- 
cording to the target PPD. 

The most general form of MCMC algo- 
rithms, the one we used in thi s work, is 
the M etropolis-Hastings algorithm (iMetropolisI 
[1951 . At iteration t, a new candidate 6* is 
sampled from a given instrumental distribution 
q{6*\ff'^), and is chosen as 6' - 9* with some 
acceptance probability depending on q, n, 0' ' 
and Q* . Distribution q has to be chosen care- 
fully, so that it can cover the whole param- 
eter space and also represent the true poste- 
rior distribution, in order to achieve a satisfac- 
tory acceptance rate. To this end, we chose a 
mixture of three distributions: a uniform and 
a Gaussian distribution with large variance al- 
low to scan the entire parameter space, while 



a Gaussian distribution with small variance is 
used for quick local refinements. 

3. A test case: a Cen A 

Since it is the closest star to the Solar System, 
a Cen A has been extensively observed. 
Effective temperature and luminosity, Tetf - 
5810 + 50 K and LIL^ = 1.522 + 0.030 
(lEggenberger et _al .' '2004!), r adius, RIR^ = 
1.2'24 + 0.003 ((Kervella et al.l l2003l) and seis- 
mic data are available and, because it is part of 
a binary, so is it s mass , M/Mq - 1.105 +0.007 
dPourbaix et al.l l2002l) . In the following we 
Umit ourselves to the non-seismic constraints. 

We used the stellar evoluti on code ASTEC 
(IChristensen-DalsgaardI Il982i 2007) to com- 
pute X(0). The estimated parameters are the 
stellar mass, M, and age, f* {Q - {M, f*}). 
By neglecting diffusion and mixing processes 
other than convection, we can fix the metallic- 
ity (Z = 0.027). The mixing-length parameter 
is assumed to be solar. 

We present in Table [T| and Fig. [T] the re- 
sults of three MCMC simulations using dif- 
ferent observational constraints and priors. In 
run 1 we retained X = |L, T^s} and uniform 
priors were used on M and f*. We used the 
same constraints in run 2 but the uniform prior 
on the mass was replaced by a Gaussian prior, 
Ti(M) oc exp(-[M/Ma - 1.105]V9.8 x lO^^). 
Run 3 aims at testing the effect of precise ra- 
dius measurements on the stellar parameters, 
we thus chose X = {L,Tsff,R] and used uni- 
form priors. 

Table [Tlgives the inferred values for M and 
f*. We present the mean value and the corre- 
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Fig. 1. Marginal distributions obtained for the mass, n(M\X) (upper row), and the age, ;r(f*|X) (lower row), 
of a Cen A. The columns correspond, from left to right, to runs 1 to 3 (see text for details). 



sponding standard deviation, derived from the 
marginal distributions 7r(M|X) and 7r(f*|X) dis- 
played in Fig. [T] The value estimated as the 
maximum a posteriori is also given. We can 
see immediately that there is a good agree- 
ment between the estimated values, which un- 
derlines the good agreement between indepen- 
dent observational constraints. 

We observe a small discrepancy between 
the mean and MAP values obtained from run 1 . 
This can be explained by the fact that con- 
vective cores start to appear in models with 
masses > 1.14 M©. Non-linear effects asso- 
ciated with convective cores will, as a gen- 
eral trend, lead to a wider range of stellar pa- 
rameters being able to reproduce the observa- 
tions (see, e.g., J0rgensen & Lindegren 2005, 
for an illustration of this phenomena), and 
hence to be accepted by the MCMC algorithm. 
Therefore the difference could be explained 
this way: by using the MAP value one selects 
only the model which reproduces best the ob- 



servations, whereas the mean value accounts 
for all accepted models. It is worth noting that 
this discrepancy becomes extremely small for 
run 3 and almost disappears for run 2, during 
which no models with convective cores were 
accepted. 
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